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Abstract 

We construct a model for tri-bimaximal lepton mixing which employs only family 
symmetries and their soft breaking; neither vacuum alignment nor super symmetry, 
extra dimensions, or non-renormalizable terms are used in our model. It is an 
extension of the Standard Model making use of the seesaw mechanism with five 
right-handed neutrino singlets. The scalar sector comprises four Higgs doublets and 
one complex gauge singlet. The horizontal symmetry of our model is based on the 
permutation group S3 of the lepton families together with the three family lepton 
numbers — united this constitutes a symmetry group A(6oo 2 ). The model makes no 
predictions for the neutrino masses. 
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1 Introduction 



Lepton mixing and non-zero neutrino masses are now established facts — for reviews and 
for the latest fits see [T]. The mixing angles in the lepton mixing matrix U have values 
quite different from those of quark mixing. The phenomenological hypothesis that 



U = [/hps = 



/ 2/VQ 1/V3 \ 

-1/y/E l/y/3 -l/y/2 
\ -1/V6 l/y/3 l/y/2 J 



[I) 



has been put forward by Harrison, Perkins and Scott (HPS) in 2002 [2]. At present, all the 
experimental data are still compatible with this simple "tri-bimaximal" mixing Ansatz. 

The hypothesis (pQ) has stimulated model building and the search for family sym- 
metries which might lead to U = Uups in a natural way. While it is not difficult to 
simultaneously obtain U e s = and maximal atmospheric-neutrino mixing [3], generating 
a solar mixing angle 9 12 = arcsin (l /y/3) is highly non-trivial and in general necessi- 
tates complicated models. In those models one often finds several scalar multiplets of the 
horizontal-symmetry group with vacuum expectation values (VEVs) aligned in a special 
way. To explain this peculiar alignment of VEVs one may have recourse to special scalar 
potentials, stabilized with the help of supersymmetry — see for instance [U El [6] — or to 
extra-dimensional models [7]. 

In two previous papers [HI E] we have enforced trimaximal mixing — which is a weaker 
hypothesis than tri-bimaximal mixing — through a model. We now show that, with very 
little extra effort, one can also achieve tri-bimaximal mixing along the same lines. In 
the model that we shall present here neither VEV alignment nor supersymmetry, non- 
renormalizable terms, or extra dimensions are required for obtaining U = £/hps- Besides 
enlarging the scalar sector of the Standard Model (SM) by several Higgs doublets and 
one gauge singlet, our model uses the seesaw mechanism [10] with more than three right- 
handed neutrino singlets, but in such a way that the additional right-handed neutrinos 
do not have Yukawa couplings to the Higgs doublets; then these additional right-handed 
neutrinos — in the present case there are two of them — can be exploited for imposing the 
desired mixing properties^] In our model lepton mixing originates solely in the Majorana 
mass matrix Mr of the right-handed neutrino singlets, and the number of independent 
Yukawa coupling constants of the Higgs doublets is an absolute minimum — only two. 

This paper is organized as follows. The model is presented in section [2j Variations 
on the symmetries of the model, and their connection to the renormalization-group evo- 
lution (RGE) of the light-neutrino mass matrix M u , are investigated in section [31 The 
conclusions are presented in section [H An appendix contains details of the computation 
of the 3x3 matrix Ai u out of the 5x5 matrix Mr. 



lr This idea had already been previously used by us for tri-bimaximal mixing, but in that case we still 
needed VEV alignment and made use of supersymmetry [5] . 
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2 The model 



2.1 Fields and symmetries 

Our model is based on the SM gauge group SU(2) x U(l). The lepton sectoJl consists 
of three left-handed SU{2) doublets D a i = {yaLi&h) (« = e,/i,r), three right-handed 
charged-lepton SU(2) singlets a R , and five right-handed SU{2) x U{1) singlet neutrinos 
VuR-, ^iR = 1,2). The scalar sector consists of one complex gauge singlet x with zero 
electric charge and four Higgs doublets <p a = (0+, 0° ) T , O = (0o > ^o) • 
The family symmetries of the model are the following: 

• Three U(l) symmetries associated with the family lepton numbers L a , 

U(l) La : D aL -> e^-£) aL , a fi -> e***a R , v aR -> e*V Q *, e [0,2tt[. (2) 

The are supposed to be sq/Wj/ broken at high energy, i.e. at the seesaw 

scale [SlUl], by dimension-three terms of the types z/^C -1 ^/,, u^C^un (C is the 
Dirac-Pauli charge-conjugation matrix). 

• The S3 permutation symmetry of the e, /1, r indices. We view this permutation 
symmetry as being generated by two non-commuting transformations: 



The cyclic transformation 



a 



< 



( D eL 

en - 



V\R 

I x- 



Hr —> t r — > e R , 

* V iiR v tR v eR , 



f/i ^ Vt ^ H^ei 
UVlR, "2R 



U 2 V2R, 



(3) 



where u = exp (2i7r/3) is the cubic root of unity with the properties cj 2 = uj* 
and 1 + + uj 2 = 0. 

The /i-r interchange [3] 



< 



D^ L <-+ D rL , 
Vr Tr> 

L X <-> X*- 



(4) 



2 We neglect the quark sector, which is immaterial for our purposes. 
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It is clear that the fields with a indices form triplet reducible representations of S3, 
while 




transform under S3 according to the complex version of the doublet irreducible 
representation, previously used for instance in [T2] H The cyclic transformation C eilT 
is softly broken by dimension-two and dimension-one terms in the scalar potential, 
but it is preserved by all the dimension-three (and, of course, dimension-four) terms 
in the Lagrangian. The symmetry J MT is not allowed to be softly broken. The VEV 
v x = (x) breaks C efMT spontaneously, but it preserves J MT because it is real; this 
is a consequence of the 7^-invariance of the scalar potential, as will be shown in 
subsection 12.31 At low energy, both C efJiT and are spontaneously broken because 
all three vacuum expectation values (VEVs) v a = (0°) o are different (see below). 

• Three Z 2 symmetries [51 [T3] 

^ ■ OR -» -OCR, 4>a -> -0a, (5) 

for a = e,/j,,T. The Z^ are supposed to be softly broken at low energy, i.e. at the 
electroweak scale, by dimension-two terms of the types 0„0/3 (a ^ Q 0o- The 
symmetry Zg i s spontaneously broken when 0° acquires the non-zero VEV v a . 



2.2 Lagrangian and lepton mixing 

The Yukawa Lagrangian has dimension four and therefore respects all the symmetries of 
the model. It is given by 

^Yukawa = ~V\ D aL a R <p a (6a) 

Ot = e,fl,T 

(6b) 

at=e,n,T 

+ f (X VmC- X vm + X* vIrC'^r) + H.c. (6c) 

The symmetries are instrumental in ensuring that only the doublet a couples to 
a.R — line fl6a|) — and that only the doublet (fi couples to the three v aR — line f l6b|) . The 
family-lepton-number symmetries ?7(1)l q are also important to enforce Yukawa couplings 
diagonal in flavour space [3] . Note that the number of Yukawa coupling constants of the 
Higgs doublets is an absolute minimum — just yi and 7/2. 

Upon spontaneous symmetry breaking (SSB) the charged leptons acquire masses m a = 
\yiv a \. Since those three masses are supposed to be all different, the scalar potential must 

3 If one wishes one may separate x m t° its real and imaginary parts, which transform under S3 
according to the real version of the doublet irreducible representation. 
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be rich enough that the VEVs v a turn out to be all different. Also upon SSB the neutrinos 
acquire, from line (I6bp . Dirac mass terms 



v eR V iiR V tR U 1R V 2R 



M 



D 



+ H.c, 



(7) 



where 



M. 



D 



{ 


a 








\ 







a 















a 

















\ 











/ 



a = y 2 v , v = <</>o) 



In the Lagrangian there are also bare neutrino Majorana mass terms. These terms 
have dimension three and are, therefore, allowed to break the family lepton numbers, but 
not the permutation symmetry S3. They are 



C 



M, 



Majorana 



+M* {vJ R C-\ R + vl R C- x v TR + v T TR G- x v eR ) 
+M 2 \v[ R C~ l [y eR + uVfj, R + u 2 v tR ) 

+M* i ul R C~ 1 u 2R + H.c. 



(9a) 
(9b) 

(9c) 
(9d) 



Together with line fl6cl) upon SSB, /^Majorana generates the neutrino Majorana mass terms 



- - (^v eR v^ R v tR u 1R v 2R ) M R C 



eR 

-T 

Kb. 

r.T 



V, 



R 



1R 

\ "Ir ) 



+ H.c, 



(10) 



where the symmetric matrix M R is 



M, 



R 



( M 
M l 
Mi 
M 2 
\ M 2 



Mi 
M 
Mi 
oj 2 M 2 
uM 2 



Mi 
Mi 
M 
uM 2 
u 2 M 2 



M 2 
u 2 M 2 
uM 2 
M N 
M 4 



M 2 
cuM 2 
uj 2 M 2 



\ 



M 4 
M' N J 



M, 



N 



* * 



M' N 



We now derive the effective light-neutrino Majorana mass terms 

\ 



2 ( uIl 



T T 

Kl Kl 



C~ X M 




H.c, 



(12) 
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where 



M, 



mIm^m 



D 



(13) 



according to the seesaw formula [10] . Because of the special form of Md in equation (J5J) , 
only the 3x3 upper-left submatrix of M^ 1 matters. One finds (for details see appendix A) 



/ x + y + t z + uo 2 y + out z + ojy + uo 2 t \ 



z + oj 2 y + ojt x + ujy + ujH 



z + y + t 



(14) 



y z + ivy + u 2 t z + y + t x + uj 2 y + ut J 
Equations (JSZl ISHj) with M 2 = M 3 tell us that 

(y,t) cc {M' N ,M N ). 



(15) 



Therefore, y/t = t> x /f*. We now make the crucial assumption that the VEV v x is real. 
This is not an unjustified assumption since it simply corresponds to the conservation of 
the symmetry 1^ by the VEV of x- It follows from this assumption that t = y, hence 



M, 



(16) 



( x + 2y z — y z — y > 

z — y x — y z + 2y 
\ z-y z + 2y x-y J 

This is precisely the M. v corresponding to tri-bimaximal mixing. Its diagonalization reads 



diag (fi 1: /j, 2 , fi 3 ) 
Hi = x + 3y — z, 
yU 2 = x + 2z, 

jj 3 = X -3y- Z . 



(17a) 
(17b) 
(17c) 
(17d) 



The light-neutrino masses are given by rrij = \fij\ (j = 1,2,3). The matrix M. v has five 
parameters, corresponding to the three neutrino masses and the two Majorana phases, 
which are completely free. 



2.3 Scalar potential 

We have demonstrated that our model leads, under the sole assumption that the VEV v x 
is real, to HPS mixing. In order to check that a real v x is viable, we proceed to analyze the 
scalar potential V of the (fi m (m = 0, e, fi, r) and x- The potential must respect both the 
three symmetries l}^ an d the permutation symmetry S3, except for the dimension-two 
and dimension-one terms, which are allowed to break softly both the Zjj, and C eflT , but 
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not I^ T . Therefore, 



V 



Ai 

+A 3 (0+0 e 0^ + 0^ <\)\<\) T + <f>l<f> T 4>l<f> e ) 

+ A 4 O O ((fricfre + 0*0,, + 0^ r ) 

+ A 5 0],0 e + 0j0r + 0r0e 0l0r) 

+ A 6 00 (0e0e + + ^T0r) 00 

»W 2 +(M 2 +(^e) 2 +H.C.' 



+A 7 
+ 



foe 



+ 



o r 



+ H.c. 

i2 



Ag (0l0e + 0^0/, + 0^0r) + Al O O 0O |X| 

+An |x| 4 + tfi (x 3 + X* 3 ) + Ail Ixl 2 + /i2 (x 2 + X* 2 ) + ^ (x + X*) 

+ Al2 [X 2 (0i0 e + ^ 2 0j,0M + ^0t0r) + X* 2 (0i0 e + ^0„ + ^0r)] 
+^2 [x (0l0e + ^0l0 M + ^ 2 0^0r) + X* (0l0e + + W<#0 T )] 

/ A*3 /ig /is /is \ / 0o \ 

//g ^ 4 /i 7 /i7 

/*8 A*7 /i5 /i6 
\ /i| /i 7 /i 6 /i 5 / 



+ 



V 0. / 



(18a 

(18b 
(18c 
(18d 
(18e 

(18f 
(18g 

(18h 

(18i 
(18j 
(18k 

(181) 



The only parameters in V which may be complex are A§ and ^7,8,9 ■ Notice the terms \i2 
and f] in line (118ip . which break C eA1T softly, and various terms in line (11811) which break 
the (and C eMT ) softly All these terms, though, preserve J^. The soft breaking of 
the I12 in li ne dUI]) is needed in order to prevent the appearance of Goldstone bosons if 
A7 = Ag = (see later). 

We want both v x and the mass of x to be at the high (seesaw) scale, while both the v m 
and the masses of the m components should be at the low (electroweak) scale. Therefore 
we must fine-tune A 12 and $2 in lines (18j ) and (I18kl) . respectively, to be extremely small, 
lest they pull the masses of the a components up to the seesaw scale@ Once A 12 and 
$2 have been tuned to be very small, the phase of v x becomes determined only by the 
terms in line (I18ip . It is clear that, if \xi is chosen negative and the product d\r] is chosen 
positive, then the minimum of V will be obtained for a real v x , with sign opposite to the 
one of $1 and rj [9J. We have thus shown that there is a range of parameters of the scalar 
potential for which the symmetry is preserved by the seesaw-scale vacuum, i.e. for 
which v x is real. 

At low scale is spontaneously broken by |t> M | 7^ \v T \. Writing 

(I^Vl ' \ v r\) (sin#, cos9) , 



This fine-tuning is a weak point of our model, but most (non-supersymmetric) models with a very 
high scale suffer from the same drawback. 
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and assuming all VEVs and coupling constants to be real, we verify that the vacuum 
potential is, as a function of 9, of the form 

a + b sin 2 20 + c sin 28 + d \A + sin 26 , 

where c cx /i 6 and d stems from the terms. Is it clear that a vacuum potential of 
this form in general leads to a non-trivial value of 6, which may moreover be very small 
if both c and d are chosen much smaller than b > 0. 



3 Variations on the symmetries and 
renormalization-group invariance 

The group structure of the model: All the symmetries of the model, and their 
respective breaking mechanisms, have been listed in section 12.11 and in principle it is not 
necessary to detail the group that they generate. Still, elucidating the group structure of 
the model may be useful for understanding the terms allowed in the Lagrangian. Following 
for instance the reasoning in [ll], the symmetry group G of our model may be described 
as the semidirect product 

G = (N x H) x S 3 , (19) 

where 

N = l}^ 1 x Z 2 ^ x Z2 is generated by the Z 2 symmetries of equation (J5J), 

H = U(l)z, e x U(l)i x U(1)l t is generated by the family lepton-number symmetries of 
equation fl2]) and 

the permutation group S3 is generated by the cyclic permutation C etlT of equation (j3J) 
and the transposition of equation 

The semi-direct product is non-trivial since neither the nor the U(l)L a commute with 
C efJiT and I^ T . The elements of G can be written as triples (n, h, s), where n G N, h G H 
and s E S3. The multiplication law of G is the usual one for semidirect products: 

(ni, hx, s x ) O2, h 2 , 8 2 ) = (n 1 s 1 n 2 s^ 1 , hxs^s^ 1 , sxs 2 ) . (20) 

In terms of 3 x 3 matrices, n is represented by a diagonal sign matrix, h is represented 
by a diagonal phase matrix and s is a matrix in the defining triplet representation of S3. 
According to section 12. 1[ the representations of G that we employ in our model are 

1 for O , 

ns for (4> e , M , r ) , 

nhs for {e R ,/i R ,T R ) , (21) 

/is for (D eL ,D^ L ,D TL ) and (z^, 1/^, ^ rii ) , 

^2(5) for and (x,X*), 
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4>a 


T 


T 




T 


T 2 



Table 1: Transformation of the multiplets under the symmetry T. The multiplets not 
shown in the table transform trivially. 

where the two-dimensional irreducible represention (irrep) of S3 is denoted D 2 (s). It is 
easy to convince oneself that all the multiplets in the list (I2T|) constitute irreps of G. 

The group G contains all the family symmetries of the dimension-four terms of the 
Lagrangian. As discussed in detail in section El there is a sequence of soft-breaking steps 
which can be described as 

G^NxS^zW, (22) 

where l^ r) is the 1 2 group generated by I„ T . 

The variations on the symmetries in the following paragraphs will only concern the 
normal subgroup H of G. 

The symmetry group A(6oo 2 ): If we remove from the three U(l) i a the global U(1)l 
associated with the total lepton number L = L e + + L T , then the normal subgroup H 
of G reduces to the set of matrices 

U (/?, 7) = diag (e 4/3 , e% e"^) , 0, 7 e [0, 2tt[ . (23) 

In this case, H x S3 is the group A(6oo 2 ), or rather a faithful irrep thereof — see [TS] for 
a study of this group. Therefore, G = N x A(6oo 2 ). 

Switching to A(54): A(54) is the group A(6r 2 ) with r = 3— for details see [13 ESI 

In this variant of our model we do not use the symmetries C/(1)l q . Instead, we define the 

matrix [8] 

T = diag (l, lu, lu 2 ) , (24) 

and use a symmetry under which the multiplets transform according to table [TJ The trans- 
formation T, together with the 3x3 permutation matrices, generates a three-dimensional 
irrep of A(54). Notice that this group is a priori smaller — hence less powerful — than 
A(6oo 2 ), but we enhance its power by allowing it to act non-trivially on the <f> a . It is easy 
to check that the Yukawa Lagrangian of equation ([6]) is invariant under T, but we still 
need the symmetries to remove from £yukawa possible non-flavour-diagonal terms [8] . 
So the symmetry group of our model is now G = N x A(54), which is finite and has 
8 x 54 = 432 elements. We may still describe G through equation (|19|) . with H replaced 
by 

H = {Amg(u p ,u q } uj~ p - q ) I p,g = 0,1,2}. (25) 

5 The latter paper uses A(54) for the construction of a lepton flavour model which is, however, totally 
different from ours. 



9 



Concerning the irreps, instead of list (1211) we now have 

1 for 0o, 

nh 2 s for (<f> e , M , (f> T ) , 

nh*s for (e R ,/j R ,T R ), (26) 

hs for [D^D^.D^) and {v eR , P^ R , v tR ) , 

D 2 {s) for {u 1R ,u 2R ) and (x,X*)- 

The breaking of T is assumed to be soft, through dimension-three and dimension- two 
terms. An important difference relative to section [2] is that T removes some of the 
dimension-four terms from the scalar potential, because it acts non-trivially on the <j) a ; 
one obtains a restricted version of equation ([TBI , viz. 

\ 7 = \s = 0. (27) 

Notice that, although we did not use the in building this variant of the model, 
eventually the Ui a turn out to be (so-called accidental) symmetries of all the dimension- 
four terms in the Lagrangian. 

Switching to A(6r 2 ) with r > 4: If the T of the previous paragraph is replaced by 
T = diag (1, cr, a*) , a = exp (2wr/r), r > 4, (28) 

then 

# = {diag (a p , <x 9 , a~ p ~ q ) | p, g = 0, . . . , r - l} (29) 

and xi S3 is isomorphic to A(6r 2 ) with r > 4. All the previous remarks, including 
tabled], still hold in this case, but there is a noteworthy exception: now we do not need to 
impose the symmetries Z^"' ', which become just accidental symmetries of all the terms in 
the Lagrangian with dimension larger than two. Eventually, the family symmetry group 
of the model is again of the form of equation (TT9]) . with G = N x A(6r 2 ) having 48r 2 
elements. 

Renormalization-group evolution of M. v : We proceed to the study of the RGE of 
M. v from the seesaw scale down to the electroweak scale. We first note that the two real 
degrees of freedom of the scalar gauge singlet \ are assumed to be heavy. Therefore, the 
renormalization-group (RG) equations relevant for the determination of Ai u at the low 
scale are simply those of a multi-Higgs-doublet SM. Those equations were derived in [18] . 
It was shown in |9] that the form of the Yukawa couplings of the charged-lepton fields — see 
line ( l6a|) — remains unchanged; only the value of y\ evolves with the energy scale. In the 
same paper [9], the importance of the quartic scalar couplings for the RGE of Ai u was 
investigated; the following sufficient conditions for RG invariance of M. v were found: 

i) The Higgs doublet O , whose VEV v is responsible for generating M. v at the seesaw 
scale, has no Yukawa couplings to the a R . In our model, the Yukawa couplings of 
the charged leptons are given by line (16"a|) at any energy scale. 
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ii) There is a symmetry, holding at the seesaw scale, which forbids dimension-five 
neutrino mass operators involving two different Higgs doublets. In our model, that 
symmetry is constituted by the three . 

iii) At the seesaw scale there is a symmetry forbidding quartic couplings of the type 



rk^k' ) (k 7^ k') in the scalar potential. In our model, this is satisfied if some 



symmetry like T leads to the condition ff27|) . 

Thus, applying the results of our previous paper [9] to the present model, we find that, 
if equation fl2~Tj) holds, then tri-bimaximal mixing holds at all energy scales in between the 
seesaw and electroweak scales. According to the preceding discussion, this is possible by 
using any of the symmetry groups A(6r 2 ) (r > 3). On the other hand, using A(6oo 2 ) 
allows both A7 and Xg to be non-vanishing, and then corrections to tri-bimaximal mixing 
from the RGE of M. v are expected. Still, it is well known that such corrections can only be 
sizable for a quasi-degenerate neutrino mass spectrum [19], an observation corroborated 
by explicit studies of multi-Higgs doublet models [18] and general considerations [20] . 

S3 versus S4: In a series of papers [21 J it has been argued that the only finite group 
capable of yielding tri-bimaximal mixing is £4, or else a larger group containing S4. We 
want to make some comments on that claim. Since S4 = A(24) [16], we can expect that a 
construction of our model in analogy to the usage of A(6r 2 ) with r > 3 is possible. This is 
indeed the case. We can place the D a i, the and the v a ^ in triplets of 5*4. Putting the 
cj) a in the reducible triplet representation of the subgroup 5*3 and adding to this scheme 
the symmetries in order to avoid non-flavour-diagonal couplings in ^Yukawa, we can 
proceed with the construction of the model just as in section [2j Actually, it is easy to 
see that this way of constructing the model amounts simply to the replacement of the 
U(l)L a by discrete lepton numbers: fermions with flavour a are multiplied by —1 instead 
of being multiplied by an arbitrary phase factor. In the language of equation ( |T9l) . in this 
case the family symmetry group is G = (N x N) xi S3 — for a complete discussion of its 
irreps see [H]. However, it appears to us that S 4 is not an adequate symmetry group 
for our model for two reasons. First, the full symmetry group, which is only effective in 
terms of dimension four in the Lagrangian, is much larger than S4 because its subgroup 

53 does not commute with the > therefore, S 4 misses an essential part of the symmetry 
structure of our model. Second, in the terms of dimension three, i.e. in /^Majorana, which 
are crucial for our model, the symmetry group is only 5*3, something that we had already 
advocated in [5]. In summary, in our model there is no compelling connection between 

54 and tri-bimaximal mixing. 

4 Conclusions 

In this paper we have proposed a model for tri-bimaximal mixing based on an extension 
of the SM with seesaw mechanism and family symmetries. The scalar sector consists of 
four Higgs doublets and one complex gauge singlet, while the fermion sector has, besides 
the SM multiplets, five right-handed neutrino singlets. The mixing matrix obtained at 
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the seesaw scale is exactly tri-bimaximal. The most straightforward version of the model 
uses as family symmetries the permutation group S3 together with three symmetries 
and family lepton numbers; the latter are softly broken at the seesaw scale. A slightly 
more complicated way to obtain the model makes use of a group A(6r 2 ) with r > 3. 
The most intricate part of the model is the stepwise soft symmetry breaking, which we 
have tried to explain carefully in section [2J Whether one uses 5*3 together with family 
lepton numbers or a group A(6r 2 ) does not make any difference, except for two terms 
of dimension four in the scalar potential. With A(6r 2 ) those two terms are forbidden 
and, as a consequence, in the one-loop renormalization-group evolution of the neutrino 
mass matrix from the seesaw scale down to the electroweak scale, that matrix retains its 
form and tri-bimaximal mixing remains exact at the electroweak scale. With S3 together 
with family lepton numbers there are the usual RGE corrections, which are quite small, 
however, whenever the neutrino mass spectrum is sufficiently non-degenerate. 

The main purpose of the model presented here is to show that in enforcing tri- 
bimaximal mixing one does not necessarily require VEV alignment, supersymmetry, non- 
renormalizable terms, or extra dimensions. As a further bonus, one can also obtain RG 
stability of HPS mixing. 

Finally, we want to stress that in our model there is decoupling of the mixing problem 
from the mass problem; the latter remains unsolved, since all lepton masses are completely 
free. 

Acknowledgements: We acknowledge support from the European Union through the 
network programme MRTN-CT-2006-035505. The work of L.L. was supported by the 
Portuguese Fundagao para a Ciencia e a Tecnologia through the project U777-Plurianual. 



12 



A Inverting Mr 

The 5x5 symmetric matrix 



M 



( 


Mo 


Mi 


Mi 


M 2 


M 3 


\ 




Mi 


M 


Mi 


w 2 M 2 


wM 3 






Mi 


Mi 


M 


cjM 2 


w 2 M 3 






M 2 


u; 2 M 2 


wM 2 


Mv 


M 4 




\ 


M 3 


wM 3 


w 2 M 3 


M 4 




1 



u = exp (2i7r/3) 



has non-zero determinant: 

det M = (M + 2Mi) {(M - Mi) 2 M N M' N - [(M - Mi) M 4 - 3M 2 M 3 ] 2 } . 



M 



-l 



Let us write 

P R 
R T Q 

where R is a 3 x 2 matrix and Q is a 2 x 2 symmetric matrix. Then, 

/ x + y + t z + u 2 y + out z + uy + u 2 t \ 

P = z + u 2 y + out x + uy + u 2 t z + y + t 

\ z + uy + u 2 t z + y + t x + u 2 y + ut J 



with 



(M 2 - M 2 ) {M N M' N - M 2 ) + (4M + 2Mi) M 2 M 3 M 4 - 3M|M 2 

detM ' 
(M 2 - M Mi) (M N M' N - M 2 ) + (M - 4M X ) M 2 M 3 M 4 - 3M 2 M 2 



x 



detM 



(M + 2Mi) m|m; 



iV 



detM 
(M + 2Mi) M 2 M N 



detM 
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